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Abstract 



Q . Truncated Levy flights are random walks in which the arbitrarily large 

steps of a Levy flight are eliminated. Since this makes the variance finite, 
"c^ I the central limit theorem applies, and as time increases the probability dis- 

"^ . tribution of the increments becomes Gaussian. Here, truncated Levy flights 

with correlated fluctuations of the variance (heteroskedasticity) are consid- 
ered. What makes these processes interesting is the fact that the crossover to 



P^ . the Gaussian regime may occur for times considerably larger than for uncor- 

Q ' related (or no) variance fluctuations. 

These processes may find direct application in the modeling of some eco- 
^ ■ nomic time series. 
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I. INTRODUCTION 

The ubiquity of the Gaussian probabihty distribution function (PDF) in physics and 
statistics is a consequence of the central hmit theorem (CLT) [Q, which states that the 
PDF of the sum of A^ independent, identically distributed (i.i.d.) stochastic variables whose 
variance is finite converges to the Gaussian PDF when N ^ oo. 

If the hypothesis of finite variance is relaxed, a generalized CLT still exists p|J^ : the PDF 
of the sum belongs to the family of Levy stable distributions, defined by the characteristic 
function (Fourier-transform) 

L{z) = exp[imz - 7|zr(l + ipz/\z\ tan(a7r/2))] (1) 

with < a < 2. For a = 2 the Gaussian distribution is recovered, while for a < 2 the PDF 
possesses power-law tails L(x) ~ Cx~^^^"^ which make the variance infinite. In this article, 
only symmetrically distributed stochastic variables are considered, for which m = (3 = 0. 

A Levy flight 0-|^ (LF) is a random walk in which the step length is chosen from the 
PDF of Eq. (1). Since Levy distributions are stable under convolution, the LF process is 
self-similar, i.e. the same Levy distribution describes increments over different time scales, 
provided the increments are appropriately rescaled. 

Levy flights appear in various physical problems @J^, in particular diffusion, fluid dy- 
namics, and polymers. However, because of their infinite variance and lack of a characteristic 
scale. Levy PDFs overestimate the probability of extreme events when used to model real 
physical systems, for which an unavoidable cutoff is always present [Q. The most direct way 
to make the variance finite is by means of truncated Levy (TL) PDFs ||^. The TL PDF is 
close to a Levy PDF for small argument, but it contains a sharp [|] or exponential cutoff 
in the tails. A truncated Levy flight (TLF) is a random walk in which the step length is 
chosen from a TL PDF. 

The TL PDF belongs to the basin of attraction of the Gaussian PDF: for large N the 
sum of A^ i.i.d. TL variables is Gaussian-distributed. For small A^ the central part of the 
PDF of the sum has a Levy shape, but the variance and higher moments are finite. For 
symmetric TL PDFs the deviation from a Gaussian may be quantified by the value of the 
normalized fourth cumulant (kurtosis, k) |]1[. This is zero for a Gaussian, and positive for a 
TL PDF. The crossover to the Gaussian regime is given by A^ ^ A'o oc /" with / the cutoff 
length. Under this condition, k becomes very small (see Eq. (20)). 

In this article, the TLF stochastic process is generalized to include a special form of non- 
linear dependence of the increments called heteroskedasticity. One can build processes with 



dependent increments such that the central part of the PDF of the sum still approximately 
behaves as in a Levy flight. What is remarkable and makes these processes interesting is 
the fact that the crossover to the Gaussian PDF may be pushed to values of A^ which are 
larger by some order of magnitudes than if the increments were independent. 

II. LEVY PDF IN ECONOMIC TIME SERIES 

Besides being interesting per se, these processes are of direct of relevance to one of the 
most noteworthy applications of LF outside physics, i.e. to the modeling of some financial 
time series §. From a physicist's perspective, the market is a very good example of a 
complex system ||9|-[TT| in which the mutual interaction and competition among a great 
number of agents (traders or speculators) with continually adapting strategies, together 
with the influence of unpredictable exogenous factors, usually produces an intricate out-of- 
equilibrium dynamics. 

In the classical equilibrium theory of economy it is assumed that "equilibrium" values 
for the prices exist, satisfying an aggregate, overall consistency condition (recalling the so 
called Nash equilibria solutions of game theory [|1^]). However, the complex dynamics of 



market prices does not seem to fit this classical picture. In particular, trading volume and 
price volatility are much higher than expected from the classical theory [0. From time to 



time, the market may display strong movements (crashes or boosts) which cannot always be 
understood in terms of rational reactions to incoming new information [Q . Instead, some 
of their features recall the physical concept of self-organized criticality []T3 . 



At last, trading volumes and variances of price increments change over time [|T6|, and may 
persist as low or high for long periods. The existence of such correlated variance fluctuations 
(heteroskedasticity) is difficult to understand in the framework of classical equilibrium the- 
ory, and few economic models can explain it (see e.g. Ref. |T^ where an equilibrium model 
capable to mimic volatility fluctuations of interest rates is developed) . 

A fundamental source of this complex dynamics may be found in the inductive, subjective 
and adaptive nature of the process leading agents to formulate the expectations which drive 
their actions P,p!8|-p^, an aspect whose fundamental features can be captured by more or 



less elaborated game or artificial-life models ^^ 

Aside the " microscopic" origin of the complexity of financial markets, another aspect of 
the problem is the phenomenological ("macroscopic") characterization of this complexity, 
in particular the study of the properties of economic time series. Given the low level of 
determinism of these series, the most fruitful description is in terms of stochastic processes. 



Martingale processes [|l| are the signature of market efficiency. In particular, random walk 
or Brownian-motion models have been used for a long time to model the increments of asset 
prices I^Dj. This might be understood by the central limit theorem if price changes resulted 
from the sum of many independent random contributions, which would seem a reasonable 
assumption. Indeed, empirical studies of financial time series have revealed gaussian behavior 
for long time scales, typically of the order of several days. However, it has been shown ^ 
that for short time scales the central part of the PDF is not Gaussian. It is well described 
by a Levy distribution, and therefore suggests an underlying LF rather than random-walk 
model. 



Non-gaussian scaling has been found in many economical or financial indeces ||8|j3T|-|37 
In financial time series, scale invariance can be characterized (i.e. the value of the self- 
similarity exponent a in Eq. (1) can be extracted) either by comparing the full PDF of price 
increments over different time scales, or by studying the time-scale dependence of some 
selected properties of the PDF. For example, the probability of return 

1 roo 

PNix = 0) = - {L{z))\ (2) 

vr JO 

which depends on the time scale A^ as N~^^", was used to extract the value a ~ 1.40 
in the high-frequency (one-minute) variations of the Standard & Poor's 500 (S&P500) in- 
dex [^,0 . An alternative quantity which may be used is the so-called Hurst exponent (see, 



e.g., Ref. H). 

As with other applications of LF, while Levy distributions describe well the central part 
of the PDF for short times, the power-law tails of these distributions are much fatter than 
observed. In particular, the observed variance is finite. All the previous remarks suggest 



the TLF as the best candidate to model these series [|36| , |38 . 

Two aspects, however, remain unexplained : one is the fact that the crossover to the 
Gaussian regime occurs at much larger times than expected from the TLF model (and from 
any model with i.i.d. increments as well). For example, the one-minute PDF Pi of the 
S&P500 index increments has kurtosis ki = 43 []8|,|38[|. If the PDF at time A^, P/v, were 
Pi convoluted A^ times with itself (with kurtosis hn = i^i/N), Gaussian behavior would be 
expected for N ^ ni (see Fig. 1). However, the central part of P/v displays Levy behavior 
up to at least A^ = 1000 [g. 

The second aspect not accounted for is heteroskedasticity: even if linear correlations are 
almost zero, there are correlations of the squared increments, i.e. the increments at different 
times are uncorrelated yet not independent random variables. Put another way, one cannot 
factorize the joint probability density of the increments at different times into the product 



of reduced densities. 

The model proposed in the following can account at the same time for both these aspects. 

III. MODELS FOR HETEROSKEDASTICITY 
A. Gaussian-type models 

Many efforts have been devoted in the past two decades to the study of time-varying 
variance, and various models have been put forward by econometricians. Roughly, two 
great classes of models exist : Auto-Regressive-Conditional- Variance (ARCH) type models 
and Stochastic- Variance (SV) type models. 

Models of both classes are usually set up in a Gaussian framework: if time is discretized 
with an elementary time step r, the increments at the k-th time-step (i.e. at t = kr) are 
assumed to be random variables of the form 

Xk= fik + Vv^^k- (3) 

Here fik is the time- varying mean of the stochastic process {fik is very small, and can safely be 
set to zero, for the short time scales considered here), Vk is a random variable representing the 
time- varying variance of the process, and the e^ are independent Gaussian random variables 
with zero mean and unit variance, independent of Vk- 

ARCH and SV models differ in the way the ak process is specified : in ARCH-type 



models |39,40] the variance at time-step A; is a deterministic function of the past squared 



increments and variances, while in SV-type models ^1| the variance is not completely de 



termined by the past data, since it contains a random contribution. With suitable choices 
of their parameters, both type of models can account for the heteroskedasticity and positive 
kustosis (leptokurtosis) of the PDF of financial series, although usually they fail to capture 
all aspects of the data |H2|. In particular, it has been shown by numerical simulations in 



Ref. IQ that the simplest ARCH-type models do not yield the Levy-type scaling of the 
PDF described above, since already at short times the value of the scaling exponent is close 
to the gaussian one. 

B. General models 

Let us assume the Xk to be zero-mean random variables with variance Vk, with a sym- 
metric PDF Pl^{xk) depending on the parameter Vk in an as yet unspecified way. The 



parameter v fluctuates in time following a stationary process. Let PAr(fi, ■■, f at) be the joint 
PDF of the variances at the different times, and let us assume the joint PDF of increments 
and variances 

i 

i.e. the increments conditional to a certain set of variances vi,--,vn are independent vari- 
ables. However, P^i^xi, ■■, xn', vi, ■■, vn) is not directly observable. The object of measurement 
is the unconditional PDF 

Pn""'{xi, ■■, xn) = I dvi- -dvNPNixi, ■■,xn;vi,--, vn), (5) 

which is only factorized if p^^Vi,--, v^) is. 

A special case of process with PDF given by Eq. (4) is 

Xk = a/^^, (6) 

where the C,k are independent random variables with zero mean and unit variance, indepen- 
dent of Vk, and with PDF -Po(0- I^ ^^is case, 

Piix) = -^Poi^). (7) 



For this special process the PDF is assumed to change with time only through a time- varying 
scale factor. The process reduces to that of Eq. (3) if Pq is Gaussian, however Pq does not 



necessarily need to be gaussian. For example, in Ref. |Q ARCH-type processes with a 
i-Student Po are used. 

The characteristic function of P^ix) is 

Pl{^z) = I dxexp{izx)Pl'{x) (8) 

^PoiV^z). (9) 

Here and in the following —>■ indicates the result for the particular case of a scale-factor-type 
process, Eq. (6). Po{z) is the characteristic function of Po{x). 

Pi{z) can be used ^ to calculate moments and cumulants of Pi{x) of any order n, 

<(i) = Hr^Pii^)u=o - v^< (10) 

Jn 

<(1) = (-if j^logPi^W|.=o -^^c°, (11) 



dz 



m^ and c^ being moments and cumulants of Pq. The assumed symmetry of P{'(x) implies 
that all its odd-order moments vanish. The variance m2{l) = 02(1) = v. The normalized 
cumulants are 

v:ii) = ^ - c°. (12) 

f 2 

The kurtosis k"(1) = ^74(1) is zero for the process of Eq. (3). 

1. Unconditional PDF for N=l 

For the Gaussian-type processes of Eq. (3), it is easy to extract piiy) from a given 
measured P^^^{x). Using Eqs. (4), (5) and (9) its characteristic function P^^'^lz) is found 
to be 

Pr{z)=Jdvp,{v)expi-'^z'). (13) 

Since pi{v) = for f < (the variance has to be positive), Eq. (13) is a Laplace transform. 
Setting z^/2 = s, P""'^(-\/2i) = C{pi(v)). Thus, the PDF of the variance pi{v) giving an 
observed P^^^{z) is the inverse Laplace transform pi{v) = C~^{Pi'^^{^/2s)). If P^^^{z) is 
a symmetric Levy PDF with index a, it follows that pi{v) is itself a Levy PDF Eq. (1) of 
index a/2 and asymmetry /? = 1, m = [§. This latter has an essential singularity at 
X = 0. If Pi"''^{z) is a symmetric TL PDF, the singular behavior close to a; = remains, but 
the decrease of pi(f) for f — ;► cxo changes from algebraic to exponential. 

For a general process such Laplace-transform relationship is lost. By differentiating 
punc^^^ at z = 0, moments m™'^(l) of P^"''^{x) of any order n can be expressed in terms of 
moments of Pi{x) : 

mr(i) = hir£^priz)u=o = «(i))., (w) 

where {f{v))v = J dvpi{v)f{v) and rn^{l) is the ra-th-order moment of Pi{x), Eq. (10). The 
kurtosis of P™'^(x) is 

;,«-(l) = ^^^"'/\^^^ - 3 ^ ^(«:o + 3) - 3, (15) 

V J (mr'=(l))2 (f)? ^ " ^ ' ^ ^ 
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where kq is the kurtosis of Pq. As is well-known, a fluctuating variance (for which {v )y 7^ 
{v)l) can produce a non-Gaussian PDF (k""''^(1) 7^ 0) even when Pq is Gaussian (kq = 0). 
Thus, a given measured P""''^(a;) may be consistent with many different choices of the couple 
{P^{x) , pi{v)), which explains the existence of many econometric models. For example, 
the observed leptokurtic character of P™^(x) may arise either from a leptokurtic Pq or from 
fluctuations of the variance or from both effects. 



2. Unconditional PDF for N > 1 

At time A^ > 1 the characteristic function of the unconditional PDF of the sum, P™^(x), 
X = J2iLi ^i, is found from Eq. (4), 

N 
Pr{z)= Uv,--dVNPN{vi,-.VN)X{P^i{z). (16) 

■^ j=l 

For independent variance fluctuations, PAr(fi, ■■, vjq) = YiiPii'^i), and P™'^(x) is simply given 
by Pi'^'^{x) convoluted N times with itself. In this case cumulants, including the variance, 
scale as A^. The kurtosis decreases as k^"'^{N) = N~^k,^"''^{1). 

If fluctuations of the variance are correlated, by differentiating Eq. (16) at z = 0, moments 
of P'^^{x) of any order n may be expressed in terms of averages of products of moments 
of Pi{x) taken at different times, the average being made over pn{vi, --jVn)- The variance 
scales as 1712"''^ (N) = Nm2"''^{l), just as for uncorrelated (or no) variance fluctuations. The 
kurtosis depends on linear correlations of the variances 



9id) 



(N) = k"""(1)/A^ + K"""(Ar) (17) 

N 

(N) = 6(2 + /?"""(!)) ^(l/AT - d/N^)g{d) 
d=i 

{ViVi+d)v - {V)l {{xlxl^^)) - ((x2))2 



where (..)t, is an average over variance fluctuations, and ((..)) is an average over fluc- 
tuations of the Xi and the variances. The normalized two-times autocorrelation of the 
squared increments, g, determines the degree of persistence of variance fluctuations. The 



simplest models of the type of Eq. (3) (such as those considered in Ref. |^) have a positive 
g{d) ~ A'exp(— (i/(io)- In Ref. |3^ g{d) ~ god'^'^'^ is found from the 5-minute increments of 
the S&P500 index. In any case, in presence of positive variance correlations, k'^"''^{N) may 
decrease with A^ much more slowly than ii g{d) = 0. Thus, roughly speaking, the slowing 
down of the decrease of k""^(A^) pushes the onset of a Gaussian regime {k^"''^{N) -^ 0) 
to much larger values of A^ than expected from independent (or no) variance fluctuations. 
The problem is that for most models, and in particular those of the type of Eq. (3), these 
heteroskedastic contributions to P^^'^{x) may be inconsistent with Levy behavior close to 
x = 0. 

To make the variance autocorrelation g explicit in P™^(z), it is convenient to character- 
ize the stochastic process followed by the variances by its multivariate characteristic function 
PN^ki, ■■, k]y). For independent Vk, Pn^^i, ■■, k]\r) = YiiPiih), Pii^) being the Fourier trans- 
form of pi(f). When correlations are present, Pn^^i, ■■, /c^v) may be expressed as 
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PNih, ■■, kn) = exp{-A^g{i-j)kikj + S)Y[pi{ki) 

i<j i 

= l[Plik^)i^ -Aj2g{z- j)hk, + ■■). (18) 

i i<j 

Here A = {m1{l))y — {v)l, and S = o{{\ki\ + ■ ■ +|fcAf|)^) contains contributions from mixed 



cumulants oi Pn{vi, ■■, v^) of order higher than two |4^, e.g. terms of type kikjki or kikjkikg 
with the time indices z, j, /, s not all equal. These describe correlations of the variances of 
higher order than those described by g. The approximation of putting S* = corresponds 
to make a Gaussian-like decoupling of these high-order correlations into products of linear 
two-times {g{d)) and equal-time correlations (e.g., {vfv2) = 3g{l){v'^) + {v){v^)). 

Let us consider, for example, models of the type of Eq. (3), which are the most commonly 
used. For these models it has been shown (Eq. (13)) that it is possible to choose pi{v) such 
that P™^{x) is a TL PDF. By using Eq. (18), g{d) may be made exphcit in P'^'^{z) , 
Eq. (16). By a few simple manipulations, it is possible to sum up all the contributions of g 
to P'^'^^z) to obtain 

Pr{z) = Pr{zfexp{—-j:g{z-j)+0{z')) (19) 

i<j 

It is peculiar to the processes of Eq. (3) that the contributions 0{z^) do not depend on g, 
but only on higher-order correlations {S in Eq. (18)). More specifically, mixed cumulants 
of order p in Eq. (18) give terms 0{z^^) in the exponent of Eq. (19). To produce a Levy 
scaling of P^^^{x) for small x, P^"''^{z) has to behave as [P""'^(z)]^ ~ exp(— A^7|2;|") for large 
z. Then, constraints should be put on variance correlations of any order to ensure that as 
z ^ oo the exponent in Eq. (19) does not diverge more strongly than z". It is also possible 
that for large z the contributions of variance fluctuations to P'^'^{z) be not determined by 
its small z behavior, i.e. by its moments [0. In turn, this may happen ii p^ivi, ■■■,vn) itself 
is not uniquely determined by its moments, i.e. PN{ki,--,kN) has singularities. Anyway, 
such non-analytic contributions would not yield a Levy scaling in general. 

IV. TRUNCATED LEVY FLIGHTS WITH HETEROSKEDASTICITY 

It has been shown above that if the elementary PDF Pi{x) is Gaussian, Levy behavior 
of P^^'^{x) may be obtained only by making ad hoc assumptions about pi{v) and the mul- 
tivariate structure of variance fluctuations, which looks somewhat artificial. It is simpler 
and more intuitive to assume, instead, that for some f , P{'(x) itself is equal or close to a 
TL PDF. There are many possible ways to choose how the parameter v should enter. If a 
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TL PDF with exponential cutoff at / ~ A ^ is used for P^, its variance and kurtosis are, 
respectively, 

^~ |cos(7ra/2)| ' " " \H ' ^^^> 

If the scaling exponent a is fixed, a time-varying v may occur either through variations of 
7 (these reflect variations of trading volumes in economic applications), or variations of the 
cutoff A or of both. 

Although the central part of a superposition of TL PDFs with varying 7 does not have an 
exact Levy behavior, this behavior still exists approximately for many choices of pi(f), i.e. 
one can find a Levy PDF with the same a and an effective 7 which describes approximately 
the central part of the superposition. For example, if f or 7 have a lognormal PDF (given 
that f , 7 > this is a simple possible candidate PDF), the central part of the superposition 
is dominated by the 7 corresponding to the peak of the lognormal (non-Gaussian region, 
high kurtosis), while the tail for 7 —i> 00 (quasi-Gaussian region, small kurtosis) mostly 
affects the tails of the superposition. In general, the smaller values of 7 determine the small 
X behavior of P^^'^{x). On the other hand, the central part of a superposition of TL PDFs 
with varying cutoff A keeps an almost exact Levy behavior, since A only affects the tails of 
the PDF. 

In the following we focus on some qualitative features of TLF with correlated fluctuations 
of the variance. We study the simplest conceivable non trivial model, which nevertheless 
contains the relevant features of this type of processes, and has the advantage that an 
accurate numerical evaluation of P™'^(x) can be made without approximations. The variance 
is assumed to fluctuate between two possible values only, Va and f?,, Va < vi,. Thus, pi{v) = 
PaS{v — Va) + (1 — Pa)S{v — ffo). If pa ^ 0.5 ouc may view pi{v) as a stylized lognormal 
PDF, with Va representing the position of the maximum and Vb representing the tail of the 
lognormal. For P^^" we use a TL PDF with high kurtosis (strongly non-Gaussian behavior). 
Pi'' should represent a TL PDF with small kurtosis (quasi-Gaussian behavior). To simplify 
the model as much as possible, and minimize the number of parameters. Pi'' is assumed to 
be Gaussian with variance Vb- Thus, P™'^ = paPi"^ + (1 — Pa)Pi'' ■ 

For the stochastic process followed by the variances, the simplest choice is a Markov 
chain |jl|. Thus, PNiyi, --^vn) is determined by assigning the probability Pa at time A^ = 1, 
and the transition probabilities paa and pbb-, i-e. the probability that ii v = Va (or v = Vb) at 
any time N, then v = Va (or v = Vb) at time A^ + 1. If the condition that the Markov chain 
be stationary is imposed, two free parameters, i.e. Pa, and pbb, are left for the chain. These 
are to be added to the 3 parameters (a, 7, A or, using Eq. (20), a, Va, Kq) which fix P,"" and 
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the single parameter (vh) which fixes P^''. Note that for this model g{d), Eq. (17), decays 
exponentially with d. 

For the purpose of illustration, we try to apply this stylized model to the S&P500 data 



of Refs. @j3y]. Four parameters are fixed by the values a = 1.4, m2"'^(l) ~ 0.016 |^ 
K^"''^{1) = 43, P^"''^{x = 0) ~ 15 measured for the 1-min. increments of the S&P500 index 
for the period 1984-1989. A fifth parameter is fixed by the decay rate do oi g{d). Although by 
the Markov chain model it is not possible to directly reproduce the measured algebraically 
decaying g{d), we do not expect this to affect the results qualitatively as long as the kurtosis 
is qualitatively well mimicked. Thus do is fixed {do ~ 126 minutes), roughly, by demanding 
K^^'^{N), Eq. (17) (see Fig. 1), to be be equal to unity at the same value A^ ~ 3700 which 
is found from Eq. (17) if one uses for g{d) the same expression g{d) ~ 0.08t/~'^'^^ found in 



Ref. [^ for the 5-min increments of the S&P500 index future for the period 1991-1995. 



Since the time window and the time step of Ref. ||3^ are not the same of Refs. ||8|j30|l, the 
value of do is only indicative. 

One parameter (pa) remains free, and is fixed at 0.9 (any value close to one gives similar 
results). 

As it may be seen in the inset of Fig. 1, the central part of P^"''^{z) is close to a Levy 
PDF. P^"''^{z) is evaluated up to A^ = 10000 by performing 4 ■ 10^ simulations of the Markov 
chain, enough to ensure full convergence of the results. These simulations yield a numerical 
estimation of the multivariate probability distribution of the variances PN{vi,--,viy), by 
which P™^(2;) is calculated through Eq. (16). P'^'^{x) is then calculated by a numerical 
Fourier transform of P^^'^{z). 

The probability of return P'^'^{x = 0) is plotted in Fig. 1. For comparison, P'^'^{x = 0) 
is plotted for uncorrelated variance fiuctuations (in this case P™'^(z) = (paPi^iz) -|- (1 — 
Pa)Pi''{z))^), and for a TLF with i.i.d. increments and again a = 1.4, m2{l) ~ 0.016, k,{1) = 
43. As expected, for the two latter models, the onset of a Gaussian regime (P™'^(x = 0) ~ 
jy-o.s-j occurs as soon as A^ ^ k""'^(1). This would not be in agreement with observations, 
since Levy scaling, P^"^(x = 0) ~ N^^^"", is observed up to at least A^ = 1000. Instead, 
for the heteroskedastic model the kurtosis decreases much more slowly (see Fig. 1) and the 
Gaussian regime occurs for much larger values of A^. An approximate Levy scaling persists 
up to A^ ~ 2000. In the inset it is shown that such approximate scaling extends to finite 
values of x as well, P]^^^{x) ~ p^c^^jy-i/a^jy-i/a g_ Gaussian scaling is estimated to 
occur only for A^ > 30000. 
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V. CONCLUSIONS 

In conclusion, TLF with correlated fluctuations of the variance (heteroskedasticity) have 
been considered. These processes may be of relevance to the modeling of some financial time 
series. An explicit numerical calculation has been made by using for the stochastic process 
followed by the variances the simplest conceivable model, i.e. a Markov chain. Parameters 
suitable to model the behavior of the S&P500 stock index have been chosen for illustration. 

The central part of the PDF of the increments during one time step, P""'^(a;), is close to 
a Levy PDF. What makes these stochastic processes interesting is the fact that Levy scaling 
of the PDF may persist for times order of magnitudes larger than for uncorrelated (or no) 
variance fluctuations. 

It has also been shown that, using the Gaussian-type models of Eq. (3), a Levy scaling 
of the PDF may be obtained only when quite ad hoc assumptions about the multivariate 
structure of variance fluctuations are made. 
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FIG. 1. Probability of return vs time N. Triangles : truncated Levy flight with Markov-chain 
heteroskedasticity (HLF) and parameters pa = 0.9, pbb = 0.9929, a = 1.4, Va = 0.00393, Ha = 500, 
Vf, = 0.1246, (see text). g{N) ~ exp(— A^/126). Circles : i.i.d. variance fluctuations, g{N > 1) = 0. 
Squares : TLF with v = 0.016, k = 43. Lower inset : kurtosis vs time N for the HLF (trian- 
gles) and for the two other non heteroskedastic processes (circles). Upper inset : scaled PDF, 
Ps{Xs) = iV^/"PjS^"'=(Xs), Xs = xiV^/" for the HLF with various values of N compared with a 
Levy PDF with a = 1.4 and 7 = 0.0037. 
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